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Abstract 



The reduction formulas for Dirac fermions are derived, using tlie exact 
^ l' solutions of free Dirac equation on de Sitter spacetime. In the framework 

, of the perturbation theory one studies the Green functions and derive the 

i-C ■ scatering amplitude in the first orders of perturbation theory. 

^ ■ 1 Introduction 

, The Dirac equation on de Sitter spacetime has been in moving or static local 

' charts leading to significant analytical solutions [2], [9], [10], [11]. The next step 

is to find solutions for the free electromagnetic field in moving local charts. The 
' problems that arises are related to the fact that curved spacetimes have specific 

, symmetries which , in general, differ from that of Minkowski spacetime. It is also 

known that form of the field equations and their solutions on curved spacetime 
are strongly dependent on the tetrad gauge and local chart in which one works 
and for that reason we don't expect to find a general solution for field equations. 
Also is important to specify that the recent astrophysical observations shows 
^ ' that the expansion of the Universe is accelerating and the mathematical model 

5^ , that could describe the far future limit of our Universe is the de Sitter model. 

Actually majority of investigations dedicated to Q.E.D on curved space- 
times don't take into considerations scattering processes. This is because in the 
present we don't have one scattering theory on curved spacetime. Our aim in 
this paper is to derive the reduction formulas for Dirac fermions and to use this 
result for developing the scattering theory on de Sitter spacetime. We will see 
that the scattering theory on de Sitter spacetime can be reproduced from that 
in Minkowski. 

We continue in this paper our work related to the developing of perturbative 
Q.E.D on de Sitter spacetime. In what follows our attention will be focused on 
Dirac field and we shall discuss only this problem. In section 3 we will derive 
the equations of fields in interaction and we use this result for constructing 
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the reduction formalism for Dirac fcrmions, in section 4 wc use this resuh to 
derive the elements of matrix for the amplitude of transition. Our conclusions is 
summarized in section 5. The results are presented in natural units h = c=l. 



2 Plane wave 

We start with the exact solutions of the free Dirac equation on de Sitter space- 
time written in [2]. Let us write the de Sitter line element [1] , 

ds^ = de - e2-*df2, (1) 

where w is the expansion factor. We know that defining a spinor field on curved 

spacetime requires one to use the tetrad fields [1] eji{x) and e''(.x), fixing the 
local frames and corresponding coframes which are labelled by the local indices 
fb, u, ... = 0, 1, 2, 3. The form of the line element allows one to chose the simple 
Cartesian gauge with the non-vanishing tetrad components: 

e9 = e-*; = ^e"-*, (2) 

so that ejx = e^e^ and have the orthonormalization properties e^ep = riixo, 

e^ep — S'^ with respect to the Minkowski metric r] = diag(l, —1,-1, —1). 
In this gauge the Dirac operator reads [2]: 

En = + ie-^'j% + ^7°. (3) 

Now let us introduce normalized helicity spinors for an arbitrary vector p by 
notation: ^a(p), 

apCxip) = 2pXUp), (4) 

with A = ±1/2 and where a are the Pauli matrices and p =| p \. For witting 
the solutions of Dirac equation on de Sitter spacetime we set: 

k=—, u± = ^±ik. (5) 
CO z 

Then the positive frequency modes of momentum p and helicity A that were 
constructed in [2] as solution of Dirac equation Euip = using the gamma 
matrices in Dirac representation (with diagonal 7°) are: 



where H^\z) is the Hankel function of first kind. 

Since the charge conjugation in a curved background is point independent 
[8], as in Minkowski case, the negative frequency modes can be obtained using 
the charge conjugation, 

u^A^) ^ Vj;Ax) = n^°(c/p,A(x))^ (7) 
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These spinors satisfy the orthonormalization relations [2] : 

/ d'x{-9Y'^UpAx)^'U^^,,^,{x) = (8) 

d^x{-gtl''U^^y,{x)^'^V^,,^,{x) = 0, 



where the integration extends on an arbitrary hypersurface t = const and 
(—(7)^/^ = e^*. They represent a complete system of solutions in the sense 
that [2]: 

/ d^p [UpAt' ^Wi^it, X') + Vp,x{t, x)Vp:xit, X')] = e-^^*6^{x - x') (9) 

The quantization can be done considering the plane wave in momentum 
representation [2]: 

il){t,x)= V'^+^ [t, f ) + V^"^ (*, x) 
= I d^'pY. x)a{p, A) + Vp^t, x)b+ip, A)] , (10) 

where the particle (a, a+) and antiparticlc (&, 6+) operators satisfy the standard 
anticommutation relations in momentum representation: 

{a(p, A), a+{p'. a')} = {bip, A), b+ip',X)} = 6^^>S'{p- p'). (11) 

From Eq.(9) follows that the equal time anticommutator takes the canonical 
form [2]: 

{i,{t,x),i,{t,x')} = e-^^*^^5\x-x'). (12) 

In any event, these are solutions of the Dirac equation and help one to write 
the Green functions in usual manner. Moreover, from the standard definition 
of the Feynman propagator: 

SfH, t,x-x')= i{0\T[il;{x)ij{x' )] |0) = e{t - t)S'^+'^ (t,t\x~x') 

-0{t' ~t)S^-\t,t ,x- x'), (13) 

in [2] was shown that: 

[Ed - m]SFit,t ,x- x') = -er^'^^S^ix - x' ). (14) 

Finally we specify that, in general, the partial anticommutator functions 
are rather complicated since iovt^t their time-dependent parts are 
complicated [2]: 

s'^^Xt,-)! ,x-x) = i{i)'^^\t,x),i)^^\t ,x')}. (15) 

With these elements we can try to develop the reduction formalism and the 
scattering theory on de Sitter backgrounds. 
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3 The reduction formalism 



We start with the interaction between spinors fields and electromagnetic field on 
de Sitter spacctimc because a theory of free fields alone has no physical content. 
The nature of physical world is revealed to observers only through the interac- 
tions between fields. In this work we use for our calculation the same formalism 
as in [3], [5]. [6]. Wc adopt the minimal coupling corresponding to the classical 
interaction of a point charge as the prescription for introducing clcctrodynamic 
couplings. Also we will note the interacting fields with tp{x) and Afi{x), where 
the hated indices indicate label the components in the local Minkowski frames. 
This fields will satisfy one system of equations that can be obtained from an 
tetrad gauge invariant action of the free Dirac field, free electromagnetic field 
and an interaction term, all minimally coupled with gravitational field: 

S[e, V, A] = J ct'x^ 1 1 [^7"£'«^ - {D^m] 

-m^^ _ If-^F"'^ - eV^T^^aV'} , (16) 

where g = \det{g^^)\, the Dirac matrices 7" satisfy {7", 7^} = 2r]°''^ and covari- 
ant derivatives in local frames. Da = P'%D^ = + r^, are expressed in terms 
of the spin connections. 

The system of equations we obtain from action (16) are: 

ij^Daip — mij} = e^^Aail) 
(Vhf)F"^) = ei^l^i^. (17) 

It is obvious from (17) that in discussing the coupling between three fields we 
are up against a nonlinear problem of vast complexity. This system of equations 
can be replaced with one system of integral equations which contain information 
about initial conditions. To do this one select two Green functions S'^{x — y) 
and D^~{x — y) corresponding to one initial condition, which help us to write 

the solutions of system (17) like follows: 

i^{x) = i^ix) - e fd^yV^S'^ix - y)^''Aa{y)^{y) 



Aa{x) = Aa{x) - e J d'^y^D^^ix - y^yW^i^iy), (18) 

where il^{x) and Aa{x) are free fields. One can verify that the first solution from 
(18) are exact solution of the first equation from (17), applying {En{x) — m) to 
ip{x), with the observation that [En{x) — m]'tp{x) = 0, and using the fact that 
the Green function must satisfy: 

[En{x) - m]S^{x -y)= -e-3-*,5^(a; - y). (19) 



4 



The Green function for Dirac field will satisfy one equation of the form (14) 
and for that reason, tpi^) is an exact solution of the system (17). 

Equation (18) offers us the possibility of constructing free fields, which are 
asymptotic equal (at t ±oo) with solutions of system (17). Now we known 
that the retarded Green functions Sii{x — y) vanishes at f ^ — oo while the 
advanced one Sa{x — y) vanishes for t — > oo. If we write first relation from (18) 
with retarded and advanced functions, 

j4 



i'ix) = i'R/Aix) - e j (fy^Sn/A{x - y)-i'' A^iyMy), (20) 

then they would define free fields, iIjr{x) that satisiy: 

hm {ilj{x)-i,ii^A{x)) = 0. (21) 

The free fields ijjji and tjjA have mass m and are equal at t ± oo with exact 
solutions of the coupled equations and represent the fields before and after 
the interaction. Now we known from Minkowski theory that the mass of the 
interacting fields may differ from that of free fields because of the connection 
between mass and energy, and the energy of electromagnetic field. In our case 
we don't know the relation between mass and energy but it is possible that this 
difference between mass of free and interacting fields to be due to the energy 
of electromagnetic field and to the coupling with gravitational field. With this 
observation the first equation from (17) can be rewritten: 

ij^Daip ~m^p = e-f^Aatp + Smij, (22) 

where the difference 6m will be calculated when one solve the system of coupled 
equations. 

Like in Minkowski case, the free fields ipn and iJja arc defined up to a normal- 
ization constant noted with y/z2- Following the same steps like in Minkowski 
case one could define in/ out fields: 

Z2tlJin/out{x) = ^{x) + e j d'^yV^Su/Aix - y)^" A&{y)'il){y) 

+5m j dS^Sn/A{x-ymy). (23) 

The in/ out free fields defined above satisfy Dirac equation, can be written with 
the help of creation and annihilation operators and satisfy conditions: 

lim {ip{x) - y/z^1pin/out{x)) = 0. (24) 
t — >^oo 

Using Eq.(22) one could write this fields as follows: 

Vz2An/out{x) = Hx) + / d'^yV^SR/A{x - y)[ED{y) - m]i;{y), (25) 
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where the Dirac operator reads Ed = i"/"Da. 

The above normahzation rule allows us to write the definition for the cre- 
ation and annihilations operators. From (10) if one uses the orthonormalization 
relations (8), one obtain: 



The creation and annihilation operators defined above satisfy the anticommu- 
tation relations (11) and from that it follows that all the properties of free fields 
will be preserved. 

Before starting our calculations we make a few remarks about the scattering 
operator. Denoting the vacuum state |0) , the one particle states for Dirac 
fermions can be written: 



If we consider two states \in, a) and \out, 13) one define the probability amplitude 
of transition from state a to state (3 as the scalar product of the two states: 
(out, /3|m, a), this is just the elements of matrix for scattering operator. This 
operator assures the stability of the vacuum state and one particle state, and in 
addition transform any out field in the equivalent in field. 

The remaining task is to construct and studies, general matrix elements 
which describe the dynamical behavior of interacting particles. We are inter- 
ested in the transitions amplitudes for interacting particles between different 
initial and final states, that is the S matrix. As in Minkowski theory we can 
construct n— particle in and out states as in the free Dirac theory by repeated 
application to the vacuum of a+(p, X)in/out and h'^{p,\)in/out- 

Let us consider the amplitude of one process in which particles from in states 
denoted with a, together with an electron l(p. A) pass in out state, in which we 
have one electron 1 (p' , A ) and particles denoted by (3. The amplitude of this 
process can be written as: 



The first term in (28) give if one use the anticommutation relations (11) the 
amplitude of one process where the electron passes from state in in out state 
without interacting with particles: 




(26) 



a"^(p. A),:„/o«t|0) = \in/out,l{p,\)) 
&^(P, A)i„/o„t|0) = \in/out,l{p,\)). 



(27) 



{out(3,l{j)' ,\ )|ma, l(p,A)) = (out/3|ao„t(p', A )|ma, l(p. A)) 

= {outp\ain{p ,\ )\ina, l{p, A)) -I- 
{out/3\{aout{p' , >^ ) - am(p', A ))\ina, l{p, A)). 



(28) 



S^^iS{p — p'){out(3\ina) . 



(29) 
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It remains to evolve the second terms, and for this we must evolve the difference: 

aout{p',>^) - ain{p',X') = J d^xe^'^'Up^x{xh°{i'out{x) - Vm(a;)). (30) 

Using (25) we obtain: 

tpoutix) - ipinix) = j d^y\f^S{x - y)[ED{v) - m]tl){y), (31) 

where S{x — y) = Sr{x — y) — Sa{x — y). Replacing (30) and (31) in (28) and 
in addition using: 

j d'xe^^^Up^xixh^'Six -y) = iUp^xiv), (32) 

we obtain: 

{outp, l{p' , A )\ina, l{p. A)) = 5^^' 5^{p — p'){otitl3\ina) 
--^ j ^Un^^,{y)[ED{v) - m]{outfmy)\^na, lip,X))d^y. (33) 

The above method can be used to any particle from in or out state. Using 
the same calculation we find the reduction formula for one positron from out 
state: 

{out (3, l(p', A )\ina, l{p, A)) = S^^'5^{p — p){outp\ina) 
+ J V^{outl3\ij{y)\ina,l{p, X))[Eniy) - m]V-, {y)d^y, (34) 

where we note En = —i-y^Da ■ The calculations of reduction formulas for elec- 
tron and positron from in state can be done using the above method obtaining 
for the transition amplitude: 

{out(}, l{p' ,X )\ina, l{p, A)) = 5xx'S^{p — p'){out(}\ina) 
^ / ^/^{outp,l{p',x')\'tl>{y)\ina)[ED{y)-rn\Uj;^x{y)d'^y, 

V2^2 J 

{outl3,l{p' , X )\ina,l{p, X}) = 5^^^' 5'^{p — p ){outf3\ina) 
j V^VpMED{y) - m]{out(3,l{p\X)\i^{y)\ina)d''y. (35) 

Now we can proceed with the reduction of the second particle. For this we 
suppose that we already done the reduction of the first electron from out state 
and we obtain (33). Now if one suppose that in the particles denoted by (3 exist 
another electron then /3 = /?' + l(p", a"), it follows to reduce this electron. The 
matrix element in which we are interested appears in (33): 

{out(3\tl){y)\ina,l{p,X)) = (out/?', l(p", a")|V' (2/)| ma, l(p, A)) 
= {outp! \ \aoutip" ,X' )'4^{y) +'0(y)aj„(p",A")]|ina, l(p, A)) 

-{outf3 \ij;{y)ain{p" ,X )|ma, l(p. A)). (36) 
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The last terms is the form {out (3 \il){y)\ina)5-y^ii 5'^[p" —p) ^ and corresponds to 
a process where electron don't interact with other particles and for that reason 
this term is not interesting for us. The first term in (36) give the amplitude 
that interest us. We will start with the evaluation of the sum using (26): 

aoutij^' ,X)i'{y) + V(y)a*«(p", a") = j e^"*C/^^^^„ (x)7° [V'o«t(a;)V'(j/) 

+i^{y)^in{x)](i'x, (37) 

then using (25) one obtain: 

ao«t(/', A")V(2/) + V'(y)am(/', A") = I ,„(a;)7° [V'(a;)V(2/)+ 

^PiyMx)] d^x+-^ l J e'^'^'Uy,^^,, (x)7° [-^(2° - x'^)S{x - z)[Ed{z) - m] 
xip{z)ilj{y) + 9{x° - z°)S{x - z)[Ed{z) - m]ilj{y)ip{z)] d^xd'^z. (38) 

Using (32) and the explicit form of Dirac operator Ed = i'J^dt + ie~'^*7'9i + 
^7° one obtain: 

aout{^',X"my)+i'{y)ain{^',X") = I e'^^U r, .„{xh'>{i,{x),^lj{y)}d^x 
+ [ e3-*c7 , {z)[Ed{z) - m] [-9{z'' - x°)^{z)^{y) + e{x^ - z^)^l^{y)^P{z)\ d^z 

V^2 J P '■^ 

+ ^ / e3-*t7 7, [z) [i7°^(^° - xO)V(^)V(y) + n°5(^° - a;°)V(2/)V(^)] d^z. (39) 

V^2 J P ' 

The last term appears because we place Dirac equation in front of parenthe- 
sis, thus acting on distributions 9 and generating the last term with changed 
sign. Integrating after dz'^ in the last term from (39) and changing the spatial 
variables z ^ x, one obtain just the first term from (39) with opposite sign. 
Now observing that in (37) the spatial integral is done at the time x^ which 
is arbitrary, one could choice = ip and obtain the chronological product 
T[xp[z)'ip{y)\. With this observation we finally obtain: 

ao«t(p", a")V'(2/) + 'ip{y)ain{p",\") = 
--^ J V^Ur,^^„{z)[Eniz) - m]T[^(^)^(j/)]rfS. (40) 

The matrix element will be: 

{outp \aout{p" A )V'(2/) +V'(2/)am(p") A )|mQ;, l(p. A)) = 



\fzi } 



^. ,/=^i77, ..(^)[^^i3W-w]W/3 |r[V'(^)V'(y)]|ma,i(p,A))(iS. (41) 

From the above calculations one sees that the reduction of the second par- 
ticle is done using the same method as for the first particle. In the matrix 
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element that interest us (41), appears two field operators that are multiplied in 
chronological order. Repeating the calculation for other particles from in and 
out states one observe that the reduction calculus is the same, indifferent what 
type of particle is reduced. 

Also we can obtain a generalization of the above formulas, supposing that 
we have n — out and m — in , Dirac particles, after we apply the reduction 
formalism finally obtain: 



The sign (— )™+" is governed by the number of sign changes dictated by the 

definition of time ordering for fermion fields. 

Now we are in the position of making important observations about the 
above reduction method. One can show that all particles, will be replaced by 
formulas that was obtain in the reduction of one particle. Also when more 
particles are reduced, in matrix element appears the time ordered products of 
corresponding field operators. Every particle will be replaced after reduction 
with expressions which depend on field operator ip{x) ((33), (34), (35)). After we 
reduce all particles from in and out states we arrive at a vacuum expectation 
value of time ordered product of fields. We don't write explicitly the spino- 
rial indices because is obviously that Eoiy) will act as (4 x 4) matrix and as 
differential operator just on spinor ip{y). One observe that thought reduction 
method the amplitudes was written as function of fundamental solutions of the 
free Dirac equation on de Sitter spacetime, and as function of vacuum expec- 
tation value of time ordered product of fields. The vacuum expectation value 
of time ordered product of fields, together with normalization constants {^J^) 
define the Green functions of the interacting fields. One can associate one Green 
function to any process of interaction. After we apply the reduction formalism, 
the Green functions of interacting fields must be calculated. 

At the end of this section we write the reduction rules for particles and 
antiparticles, denoting one electron by l{p, A) and one positron by l(p. A), after 
reduction of particles from in and out states one obtain: 




x(0|T[^(a;i)...^(a;„)V;(%)...V^(2/m)]|0)[£D(%) - m]U^-.^^.{yj). (42) 



l(p', A ) out 






(43) 
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4 The perturbation theory 



The Green functions of the interacting fields can't be calculated exact and for 
that reason wc will use perturbation methods. The form of the amplitudes 
obtained from reduction formalism allows one to use perturbation calculus. 

It is clear now that the entirely perturbation theory on de Sitter spacetime 
can be reproduced from Minkowski theory [3] , [4] , [5] , [6] . For calculating the 
Green functions we must write then as functions of free fields, because we know 
their form and properties. We write the Green function in generally as follows: 

G(yi,2/2,...,yn) = — i-(0|T[^(yi)V'(y2)...V'(yn),.5]|0), (44) 
{0\o\0) 

where 5 is a unitary operator and have a closer form with the same operator from 
Minkowski theory. Like in Minkowski case the operator S must be correlated 
with scattering operator S. One can show that this two operators are equal up 
to a phase factor. _ 
Then entire perturbation calculus is based on development of operator S: 

g _ J-g-' / V^Ci(x)<fx _ ]^ _|_ 

£ ^ / i-9r'^n'^Ii^l)■■■'^I{^n)]d'x^...d^Xn^ (45) 

n=l 

where the density lagrangian of interaction reads: Ci{x) = — e : 'ip{x)j°'A^{x)tp{x) 
—5m : r{j{x)ijj{x) :. Each term from (45) corresponds to a rang from perturba- 
tion theory. Replacing (45) in the expression of Green function (44) one obtain 
perturbation series which allows one to calculate the amplitude in any order. 
The term of rang n of this development is: 

-TtJtT /(-ff)"/'(0|T[V'(yi)V'(2/2)...^(yn),/:/(a:i).-/:/(a:n)]|0)rfV...d^a;„. 
n!(0|6|0) J 

(46) 

The evaluation of the integrant from (46) is the same as in Minkowski case, we 

have a cinematic part which is obtained from reduction formalism and a dynamic 
one represented by operator S. Following the same steps as in Minkowski theory 
one will make the T contractions between cinematic and dynamic part, with the 
observation that the T contractions between fields from cinematic part will not 
give contributions to the scattering amplitude. Also in the case of T contractions 
between fields from dynamic part we have two possibilities. One is that all 
fields from dynamic part coupled fields from cinematic one and the second is 
that one part of the fields from dynamic part contract between them. The 
second possibility will give one term of the form (0|6'|0), which will simplify the 
nominator. 

As an application to our formalism we can obtain the scattering ampli- 
tudes in first orders in perturbation theory, thus completing the framework that 
one needs for calculating scattering processes in the first order of perturbation 
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theory, on dc Sitter expanding universe. Now using the reduction forrnahsm 
developed in section 3 for one amphtude of the form {out, l{p' ,X )|m, l{p, A)) , 
one obtain the development: 

{out,l{p',X )\in,l{p,X)) = S^y6^{p-p') - 
j i-9)U;,^y{y)[ED{y) - m]{Om{ym^)m[ED{z) - m]U^^,x{z)d%d^z. (47) 

Using (46) and neglecting the first term in (47), one obtain the for scattering 
amplitude in first order of perturbation theory: 

Ai^^ = -^J I li-gf'U,^^,{y)[Eniy)~m]{0\T[:i,iy)^iz): 

: ^{x)jfiAi'{x)i;{x) :]|0)[^d(2) - m]Up,x{z)d''yd^zd^x. (48) 

After we make the T contractions between cinematic and dynamic part using 
the method from Minkowski theory, and use (13) and (14), finally obtain: 

Ai^f = -ie J ^/^U;, y {x)jf,Al'ix)Up-,xix)d^x. (49) 

The above expression is just the scattering amplitude that we use in our previous 
work [7], for calculate the Coulomb scattering on de Sitter expanding universe. 
Following the same steps as above one can obtain the scattering amplitude in 
superior orders in perturbation theory. With the above calculations we establish 
the general rules of calculation which can be used in the language of Feynman 
graphs. 



5 Conclusion 

In this paper we investigate the reduction formalism for solution of the free Dirac 
equation on de Sitter background. We obtain that the reduction formalism can 
be calculated using the same method as in Minkowski case. Also we show 
that after reduction of particles from in and out states one obtain the vacuum 
average of fields written in chronological order. As in the Minkowski case to 
any interaction one can associate one Green function, that can be evolved using 
perturbation theory. From our formalism of reduction and using perturbation 
theory, we deduce the correct definition for the scattering amplitudes. 

From our point of view is important to do the same calculations for the 
electromagnetic field, thus completing this theory and the framework that one 
needs for developing perturbative Q.E.D on de Sitter spacetime. 
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